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Abstract 



The Planck mission, originally devised for cosmological studies, offers the op- 
portunity to observe Solar System objects at millimetric and submillimetric wave- 
lengths. We concentrate in this paper on the asteroids of the Main Belt, a large class 
of minor bodies in the Solar System. At present more that 40000 of these asteroids 
have been discovered and their detection rate is rapidly increasing. We intend to 
estimate the number of asteroids that can can be detected during the mission and to 
evaluate the strength of their signal. We have rescaled the instrument sensitivities, 
calculated by the LFI and HFI teams for sources fixed in the sky, introducing some 
degradation factors to properly account for moving objects. In this way a detection 
threshold is derived for asteroidal detection that is related to the diameter of the 
asteroid and its geocentric distance. We have developed a numerical code that mod- 
els the detection of asteroids in the LFI and HFI channels during the mission. This 
code perfoprms a detailed integration of the orbits of the asteroids in the timespan 
of the mission and identifies those bodies that fall in the beams of Planck and their 
signal stenght. According to our simulations, a total of 397 objects will be observed 
by Planck and an asteroidal body will be detected in some beam in 30% of the 
total sky scan-circles. A significant fraction (in the range from ~ 50 to 100 objects) 
of the 397 asteroids will be observed with a high S/N ratio. 

Flux measurements of a large sample of asteroids in the submillimeter and mil- 
limeter range are relevant since they allow to analyze the thermal emission and its 
relation to the surface and regolith properties. Furthermore, it will be possible to 
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check on a wider base the two standard thermal models, based on a nonrotating or 
rapidly rotating sphere. 

Our method can also be used to separate Solar System sources from cosmological 
sources in the survey. This work is based on Planck LFI activities. 

Key words: Astronomical and space-research instrumentation, Astronomical 
observations: Radio, microwave, and submillimeter, Solar System: Asteroids 
PACS: 95.55.n, 95.85.Bh, 96.30.Ys 
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1 Introduction 



The Planck ESA mission 2 will perform a high-angular resolution mapping 
of the microwave sky over a wide range of frequencies, from 30 to 900 GHz. 
These data will have important implications on different fields, from cosmol- 
ogy and fundamental physics to Galactic and extragalactic astrophysics. The 
Solar System astronomy will also benefit from the Planck mission since it will 
offer the opportunity to perform a survey of Solar System objects at milli- 
metric wavelengths. Planck will be sensitive to the millimetric emissions from 
planets and from a significant fraction of the asteroidal population. In this 
paper we focus on the detection of asteroid radio emissions and on the identi- 
fication of the asteroidal targets that will be observed during the mission. The 
relevant parameters for the possible detection of a minor body by Planck are 
the geocentric distance and the diameter. We will concentrate in this paper on 
Main Belt asteroids that can be more easily detected since their orbits are at 
relatively small geocentric distances and several objects have diameter larger 
than 50 km. 

Radio observations of asteroids from the Earth were performed by Redman et 
al. (1998) at seven different frequencies ranging from 150 to 870 GHz. Redman 
focused on the information that radio data can give to the thermal models of 
asteroids. Their set of targets included 5 asteroids already observed by Al- 
tenhoff et al. (1994) and two new objects. Altenhoff et al. (1994) observed, at 
2 http://astro.estec.esa.nl/SA-general/Projects/Planck/ 
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250 GHz, 15 among the largest asteroids, but only to determine the diameter 

and to compare their measurements with other methods. 

The emissivity obtained at radio frequencies may probe the temperature of 

the body at different depths on the surface allowing to derive precise values 

of the average temperature. Moreover, relevant physical data concerning the 

density of the regolith layer covering the surface of the asteroids can also be 

determined. 

Radio observations of asteroids by Planck offer a unique opportunity to 
improve our knowledge on the thermal emission in the millimetric and sub- 
millimetric domain, increasing our knowledge of the physical nature of the 
surface layers of the objects. Considering that radio data provide an estimate 
of the density on a scale related to the wavelength of the observations, the ad- 
vantage of using Planck will be in the possibility of gathering data on a large 
frequency spectrum and then to probe the temperature and regolith density 
at different depths below the surface. 

Furthermore, the Planck data can help to refine the thermophysical models 
using also physical studies of the thermal properties of stony and FeNi me- 
teorites as well as inferences about surface physical properties derived from 
asteroid radar studies. It will be possible to relate the Planck data to the 
two standard models which are commonly invoked to predict the thermal ra- 
diation from asteroids, considering the asteroid as a nonrotating sphere or a 
rapidly rotating sphere. 

For instance Vesta has been observed on a wide range of wavelengths, from 
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submillimeter to centimeter, and the ratio of the observed flux, at each wave- 
length (from 12/xm to 6cm), divided by the flux expected from a blackbody 
at the temperature of a nonrotating asteroid shows a behavior different from 
the standard models (Redman et al. 1992). A successful model for the ther- 
mal emission must explain why the apparent flux remains below the rapidly 
rotating sphere value, and much below the nonrotating sphere value, over the 
entire submillimeter range. Several mechanisms can lower the apparent tem- 
perature at radio wavelengths as reflections at the outer surface of the regolith, 
reduction in the actual temperature of the deeper layers of the regolith due to 
diurnal temperature variations at the surface, or scattering by grains within 
the regolith which can reduce the emissivity in a wavelength dependent fash- 
ion. The actual emission from an asteroid may combine all three effects in 
differing amounts at each wavelength, but only very few objects have been 
observed in the submillimeter and millimeter range and it is not possible to 
identify a typical behavior or to relate the results to different asteroid taxo- 
nomic classes. 

It is worth to point out that the main belt asteroids may represent a source 
of background for the PLANCK mission and, as a consequence, they should 
be considered in the definition of the scientific and technical aspects of the 
mission. 

The thermal emission of asteroids, apart from their intrinsic scientific rele- 
vance, has also to be considered as a potential source of spurious detections 
in the analysis of Galactic or extragalactic sources. The precise determination 
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at each epoch of the eventual presence of an asteroid, with its point emission, 
in some of the Planck beams is a critical requirement to avoid systematic 
errors in mapping the astrophysical sources of radiation. The high nominal 
accuracy of the PLANCK instruments suggests a tight control even for the 
contamination from asteroids below the detection threshold. 
Section 2 is devoted to the discussion of the Planck sensitivity for the de- 
tection of moving sources at different frequency channels. The averaged final 
Planck sensitivities will be rescaled in order to take into account the po- 
sitions of the different antenna beams in the telescope field of view and the 
main properties of the Planck scanning strategy. These results are used in 
Section 3 to determine the threshold to detect asteroids in the Planck data 
streams, according to their typical temperature and size. In Section 4, we 
describe our numerical code to compute the asteroid transits on the differ- 
ent Planck beams and to estimate their signals. The results of this analysis 
are presented in Section 5. Finally, we discuss the results and draw our main 
conclusions in Section 6. 



2 Planck performances 

The Planck surveyor will observe the sky at nine frequencies with two instru- 
ments having different angular resolutions and sensitivities: the Low Frequency 
Instrument (LFI) that covers the range 30-100 GHz, and the High Frequency 



Table 1 

Instrumental performances and confusion noise estimates for Planck LFI. The 
bandwidths are 20% for the LFI channels and 25% for the HFI ones. 

^eflf N ra d-bol FWHM Cr no isc Cnoisc ^CMB fcx.sou. C^Gal Cnoise,MAP 

(GHz) (unpol;pol) (arcmin) (/iK) (mJy) (mJy) (mJy) (mJy) (mJy) 
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Instrument (HFI) that covers the range 143-857 GHz. In Table 1 we report 
the relevant parameters for both LFI and HFI (see columns 1, 4 and 5). The 
relative bandwidths of the LFI and HFI are respectively ~ 20% and ~ 25% 
(see column 2 of Table 1). 

In this work we will refer to the nominal Planck sensitivity per resolution 
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element (a squared pixel with side equal to the Full Width at Half Maxi- 
mum (FWHM) of the corresponding beam), as recently revised by the LFI 
Consortium (Planck Low Frequency Instrument, Instrument Science Verifi- 
cation Review, October 1999, LFI Design Report, private reference; see also 
the LFI proposal, Mandolesi et al. (1998)) and as reported in the HFI pro- 
posal (Puget et al. 1998). The sensitivity is related to the measurement of 
microwave anisotropy in terms of antenna temperature and flux 
The normalized sensitivity per pixel in the sky (i.e., the sensitivity divided by 
the averaged sensitivity of the observed pixels in the sky) for a typical Planck 
beam improves with the module of the ecliptic latitude up to a ring of opti- 
mal sensitivity which contains a small unobserved sky region. The normalized 
sensitivity averaged over the ecliptic longitudes as a function of the ecliptic 
colatitude is close to unit at the ecliptic colatitude 9 e ~ 50° and 130° and it 



is quite well approximated by the law ysin^g/sinSO at least for ecliptic lati- 
tudes between about —70° and +70°, relevant for Solar System objects. The 
detailed behavior of the sensitivity on 9 e is clearly dependent on the choice of 
the scan angle a, the adopted scanning strategy and the beam position on the 
telescope field of view. 

The Solar System bodies are moving objects and most likely they have a 
variable flux depending from their heliocentric and geocentric distances. It 
means that their observation and estimate of their fluxes can not take easily 
advantage from the twice coverage of the sky, during the planned period of 
observation of 14 months, nor from the coadding of data from different feeds 
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if they look at different sky positions at the same time. Further details on the 
Planck sensitivity relevant for moving bodies and variable sources can be 
found in (Burigana 2000) and (Burigana 2001). For the LFI sensitivities we 
take into account the most recent computations of the LFI beam positions of 
the telescope field of view (Sandri & Villa 2001). The HFI channels can play 
a crucial role for the study of the moving bodies of the Solar System, as their 
spectra are typically close to blackbodies, or modified blackbodies, at temper- 
atures significantly higher than the cosmic microwave background (CMB) or 
are dominated by dust emission components or by molecular emission lines, 
particularly relevant at the highest Planck frequencies (see, e.g., Neufeld et 
al. (2000)). 

The Planck receivers are affected by the l/f a (with a in the range from 
~ 1 to ~ 2) noise that introduces drifts in the time ordered data identified 
as stripes in the final maps (Janssen et al. 1996). The time ordered data have 
to be cleaned from these drifts by using appropriate data reduction methods. 
Destriping algorithms (see, e.g., Delabrouille (1998), Maino et al. (2001), and 
references therein) are able to efficiently remove 1/ f a noise drifts by working 
directly in the time domain, fully preserving the time ordering of the data and 
allowing to exploit the cleaned data to properly reveal the signatures from 
moving and variable sources. Then it is possible to neglect the drift contam- 
ination in the following discussion. We note also that destriping algorithms 
(as well as other classes of reduction methods) use the comparison between 
repeated observations at different times of the same sky positions to remove 
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drifts. Therefore, cleaning the time ordered data from astrophysical moving or 
variable sources by properly flagging the time ordered data may be useful to 
improve the destriping codes, according to the S/N ratios of the contaminat- 
ing sources (the possibility to use flagged data have been already implemented 
in the most recent version of the LFI destriping code, see Stanghellini et al. 
(2001)). Running this task represents a minor effort for the data reduction 
pipeline and it can be implemented starting from the software developed with 
this work. 

In addition to the instrumental noise, the sky itself introduces a confusion 
noise given by the Galactic and extragalactic temperature fluctuations and 
from the CMB anisotropies, relevant at different angular scales, 9 ~ 180°/£, 
where £ is the multipole of the harmonic expansion of the temperature fluctua- 
tion pattern. We report in Table 1 the fluctuation level of the CMB anisotropy 
assuming a rms thermodynamic temperature fluctuation of about 95 /xK as 
derived for a standard cosmological flat model, approximately COBE/DMR 
normalized with cosmological parameters compatible with the present con- 
straints on CMB fluctuations at moderate multipoles, as derived from recent 
balloon experiments. Of course, the accurate determination of the CMB con- 
fusion noise at small scales is the Planck main goal. 

Then we report the extragalactic source confusion noise (by taking into ac- 
count both radio and infrared sources), as evaluated by Toffolatti et al. (1998) 
and revised by De Zotti et al. (1999a) and De Zotti et al. (1999b). We estimate 
a more precise source confusion noise as a ex , sou , = ^/C//(FWHM/rad) where 
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the angular power spectrum Ci is the sum of the Ci of radio and far-IR source 
fluctuations, as quoted by De Zotti et al. (1999b) for a conservative detection 
limit of 1 Jy (see their Table 6). These estimates may be then pessimistic. On 
the other hand, it would be probably difficult to subtract sources in real time 
with high accuracy; in addition, the source clustering, neglected in this work, 
may significantly increase the source confusion noise, at least for some cosmo- 
logical models, particularly at the highest Planck frequencies (Magliocchetti 
et al. 2001). 

The Galactic fluctuation levels reported here are only indicative and refer to 
moderate and high Galactic latitudes, as large variations are present in the 
sky. 

The relevant global sensitivity for point source detection/observation is typ- 
ically assumed to be the sum in quadrature of all the sources of confusion 
noise, multiplied by a proper constant (typically in the range from ~ 2 to 5). 
We report also in the last column of Table 1 the sensitivity of the MAP exper- 
iment 3 (~ 35yuK, in terms of thermodynamic temperature, on squared pixels 
with 0.3° side at each frequency channel, 22, 30, 40, 60 and 90 GHz), rescaled 
at the corresponding Planck beam size. In fact, we can argue that after the 
MAP mission the sum of the CMB and foreground temperature fluctuations 
on each sky pixel will be known with the MAP accuracy. 
CMB anisotropy experiments provide time ordered data and then maps of 

3 http://map.gsfc.nasa.gov/html/ 
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temperature fluctuations, ST, from the average temperature; the temperature 
fluctuations can be expressed in terms of fluctuations of antenna temperature, 
ST ant , or in terms of fluctuations of thermodynamic temperature, STtherm- 
Their ratio is given by: 

ST ant = x 2 exp(x) 

ST therm ~ [exp(x) - l] 2 ' U 

where x = hv/kTo, To ~ 2.725 K being the CMB monopole thermodynamic 
temperature as established by Mather et al. (1999). The same formula holds 
for ratio between the rms noise of the observed antenna temperature fluctua- 
tion, AST ant , and the rms noise of the observed thermodynamic temperature 
fluctuation, AST t h er m- 

In the study of discrete sources it is preferable to express the signal fluctuations 
in terms of flux fluctuations. For a point source with flux F v observed with 
a beam response J normalized to the beam maximum response, the observed 
antenna temperature is given in general by: 

T v , ant , obs = 10- 41 (c 2 /2A;)[(F ! ,/Jy)/(z//GHz) 2 ](J/ J JdQ) . (2) 

4-7T 

Assuming a beam quite well approximated by a Gaussian shape, as in the case 
of Planck beams, the integral of the antenna pattern is given by J 47r Jdfl ~ 
(7r/41n2)(FWHM/rad) 2 . In this work we deal with sources with weak or mod- 
erate flux, i.e. essentially detectable only when they fall within a pixel with a 
side ~ FWHM about the beam center. For a first order evaluation of the source 
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flux sensitivity estimate we can then neglect the information contained in the 
time ordered data outside the pixel with a side ~ FWHM about the beam 
center, assume the Planck sensitivity per pixel of side ~ FWHM, and approx- 
imate the beam response as unit within a pixel with a side ~ FWHM about the 
beam center and null outside in the relationship between the antenna tempera- 
ture and the flux reported in the following, i.e. J/ J 4n JdVt ~ l/(FWHM/rad) 2 
within a pixel with a side ~ FWHM about the beam center and J / J 47r JdQ = 
outside. 

In this approximation, the relationship between the flux fluctuation and the 
antenna temperature fluctuation induced by a discrete source on a squared 
pixel with side A6 ( ~ FWHM) is: 

6Bu/Jj ~ 30.7[z//GHz] 2 [A#/rad] 2 , (3) 



[5T ant /mK] 



and the same formula holds for ratio between the rms noise of the observed 
flux fluctuation, A5B V , and the rms noise of the observed antenna tempera- 
ture fluctuation, AST ant , on a pixel of side ~ FWHM. This approximation is 
used to translate the HFI sensitivities quoted by Puget et al. (1998) in terms 
of total flux fluctuation to the antenna temperature fluctuation sensitivities 
reported in Table 1. More detailed analyses, able to take into account the 
complexity of beam shapes, the pointing uncertainty and other relevant sys- 
tematical effects will be considered in future technical works and included in 
the data analysis pipeline. 
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We will derive the Planck sensitivity at the different channels, relevant for 
the study of moving objects from the Planck time ordered data, by rescal- 
ing the sensitivities of the final Planck maps. At v & 100 GHz it would be 
possible to exploit the microwave maps derived from MAP data and cleaned 
at the level of MAP sensitivity from the effect of source flux variations and 
moving objects. In this case the la sensitivity levels, relevant for the follow- 
ing discussions considering the MAP data, will be the sum in quadrature of 
the Planck sensitivities, appropriately rescaled, and of the sensitivities of 
the final MAP data. These values will be also the sensitivities to be used to 
detect and to study the moving objects before producing the Planck maps. 
On the other hand, after the Planck data analysis and the production of the 
Planck cleaned maps at all frequencies, we may be able to take advantage 
from the knowledge of the sky fluctuation in the positions of the considered 
moving objects. The sensitivity will be close to the final one of the Planck 
maps multiplied by y/2, i.e. about a factor 2 better than the MAP final sensi- 
tivity in the common frequency range, because a moving object will not have 
the same position in the following Planck sky coverage. 

2.1 Planck sensitivity for moving sources 

The Planck sensitivity relevant for the study of moving sources has to be 
properly rescaled with respect to that reported in Table 1 to include several 
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factors of sensitivity degradation. 

A first sensitivity degradation factor, ~ y/2, derives by considering only a sin- 
gle sky coverage and provides a lower limit to the sensitivity utilized in this 
work (we neglect in the present estimates the small increase of this degrada- 
tion factor with respect to the value y/2 introduced by the limited sky regions 
observed during the mission). The estimate of the Planck best sensitivity 
for moving sources assumes to be possible to properly use and to average the 
information from all the receivers of a specific frequency channel. Of course, 
a realistic situation may be less favorable, because of the motion of the con- 
sidered object. 

The LFI beams for each channel are located in a ring with a radius of about 
4°, on the Planck telescope field of view, around the line of sight. The HFI 
beams are located closer to the center and they may be also at few degrees 
from the telescope line of sight. On the Planck telescope field of view, two 
coordinates U,V are typically used to identify the beam positions. The direc- 
tion in the sky of the coordinate V is parallel to the scan circle of the Planck 
telescope line of sight, generated by the spacecraft rotation around its spin 
axis. The direction in the sky of the coordinate U, orthogonal to V, is, at least 
for simple scanning strategies (i.e., with the spacecraft spin axis always par- 
allel to the antisolar direction), parallel to the spin axis shift direction during 
the mission. 

Given the spread of the Planck LFI beams on the telescope field of view, 
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another degradation factor, ^N TiCySky) where N ric ^ sky , ranging from 2 to 16, is 
the number of radiometers per channel pointing at the same sky position at 
the same time, has to be taken into account. In column 3 of Table 1 we report 
the number of LFI radiometers (all intrinsically sensitive to the polarization) 
and of the unpolarized and polarized HFI bolometers. In the case of the LFI 
channels the total temperature measurement is given by adding the signals 
from the two radiometers coupled to a specific feed horn (or beam); therefore, 
N r ic,sky = N rad /2 = Nf eed . In all the cases, except for one feed at 44 GHz, a 
pair of LFI feeds look at the same scan circle in the sky, being located in the 
focal plane unit to follow the sky scan direction as the spacecraft spins (each 
LFI beam, except one at 44 GHz, is symmetrically located to another LFI 
beam with respect to the axis V=0, i.e. they have the same U coordinate). 
The factor, ^N ric ^ sky , reported in Table 2 at 44 GHz, provides a pessimistic 
value for our sensitivity estimates. A different factor, for most of the cases, 
and reported in Table 2 for the 30, 70 and 100 GHz channels, is \J N riCySky /2. 
Puget et al. (1998) obtained the sensitivity on temperature fluctuation mea- 
surements by combining both unpolarized and polarized detectors for the HFI 
channels at 143 and 217 GHz. At 545 GHz only polarized detectors will oper- 
ate, whereas at 100, 353 and 857 GHz only unpolarized detectors are planned. 
The HFI Focal Plane Unit arrangement is more complex than for LFI and 
it will be re-defined in few months (B. Maffei, private communication at the 
Planck Systematic Effect Working Group meeting, ESTEC, June 2001). To 
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Table 2 

Instrumental noise estimates for moving sources. The values reported in the columns 
2 and 3 refer to the Planck sensitivity for moving sources at middle ecliptic lat- 
itudes. In the case of low ecliptic latitude objects, we report also the la overall 
confusion noise including CMB and astrophysical fluctuations, neglecting the in- 
formation from future space missions (column 4), by considering the information 
provided by the MAP survey (column 5), and, finally, the sensitivity appropriate to 
a Planck single sky coverage (see also the text). 
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translate the sky map temperature sensitivity of Table 1 to the simple esti- 
mate of upper limits for moving sources quoted in Table 2 we consider that 
analogously to the case of the LFI, the HFI feeds at each frequency are aligned 

18 



in two or three groups. The feeds of each group, looking at the same scan circle 
in the sky, have a very similar value of U. In the case of the Planck HFI, 
we have to include a further degradation factor in the range from ~ a/2 to 
~ \/3- In the sensitivity estimates reported in Table 2 we conservatively use 
the value \/3. 

A last degradation factor ~ ^FWHM/A0 s , where A9 S (2.5 arcmin) is the 
spin axis shift between two consecutive hours, takes into account the possi- 
bility that a moving object falls out of the beam after a spin axis repointing. 
This factor has to be applied both to LFI and HFI. 

Taking into account all these sensitivity degradation factors, the nominal 
Planck sensitivities can be resumed in the ranges reported in column 2 and 
3 of Table 2; we provide these sensitivities in terms of antenna temperature 
and flux. 

In practice, we have to keep in mind that the interesting moving objects are 
generally located close to the ecliptic plane, where the Planck sensitivity is 
~ 1.15 worst than the averaged sensitivity. Therefore, the values reported in 
columns 2 and 3 of Table 2 have to be typically multiplied by ~ 1.15. By tak- 
ing this factor into account, we report also in Table 2 the sum in quadrature 
of the Planck sensitivity for moving objects as quoted in column 3 and of 
the astrophysical and cosmological confusion noise sources of Table 1 (column 
4), the sum in quadrature of the Planck sensitivity for moving objects as 
quoted in column 3 and of the sensitivity of MAP final maps (column 5) and, 
finally, the sum in quadrature of the Planck sensitivity to moving objects as 
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quoted in column 3 and of the Planck sensitivity appropriate to half mission 
(column 6). These numbers are indicative because of the large variations of 
the Galaxy intensity in the sky and the uncertainty of extragalactic source 
fluctuations at v £ 143 GHz (intermediate values are adopted in this work). 
An important point is represented by the opportunity to combine the millimet- 
ric observations of asteroids performed by Planck with other data obtained 
in the same period and in different spectral regions, using ground-based tele- 
scopes. 

The Planck asteroid observation at different frequencies will be obtained 
within ~ few days, because of the angular size of Planck receiver array and 
the spin axis repointing of 1° per day. Then we may predict to observe an 
asteroid with Planck within a week. 

Planck will be injected in its Lissajous orbit around L2 after about 3-4 
months from the launch and the PLANCK nominal scanning strategy will 
be certainly decided at least many months before the launch. The effective 
PLANCK scanning strategy may be slightly different from the nominal one and 
will be accurately reconstructed in few days from the telemetry data. These 
requirements will allow to organize an observing campaign, of the objects de- 
tected by Planck, on ground-based telescopes several months in advance. 
The first rough estimate of the flux for an asteroid observed by Planck will 
be delivered in a few days, waiting for a longer period for a more accurate 
value, at least 6 months corresponding to the survey duration. This estimate 
should allow a first comparison with the ground-based observations. 
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3 Threshold detection for asteroids 



In the radio range the mean equilibrium temperature T eq of an asteroid, av- 
eraged over its extent, is given by: 

T eq = f(l-p) 1/4 r- 1 / 2 , (4) 

where p is the bolometric albedo and r the heliocentric distance. In the non- 
rotating case / is 329, and for our fast rotating asteroids (few hours) / is 
277. For instance in the case of 1 Ceres, the mean distance from the Sun is 
2.767 AU and the albedo is 0.06 (Morrison 1974), which give an equilibrium 
temperature of 165 K. 

The asteroid equilibrium temperature is related to the asteroid antenna tem- 
perature, both averaged over the asteroid extent, by: 

^u,ant,obj = e v~^ eq j (5) 

with e u the emissivity of the asteroid at the frequency v. In this work we 
need to define a mean value of T* t 6 ■ for the main belt asteroids in order 
to determine the threshold detection. We will assume the emissivity close to 
unity, since we are considering objects with T eq rs 100 K and working at 
millimetric wavelengths. 

The antenna temperature observed at a specific frequency v is given by the 
convolution, over the whole solid angle, of the antenna temperature, of the 

21 



considered asteroid, with the antenna beam pattern response (we assume the 
antenna responses normalized to the maximum, i.e. J = 1 at the beam center 
direction): 

rp /,x _ / 4 Aant, 6j(7> t ) J (^^) dn 

in this equation 7 is a unit vector identifying a generic direction in the sky. 
As evident from this equation, the observed antenna temperature depends on 
the considered observation time because of the spacecraft rotation and scan- 
ning strategy and the variation of the asteroid properties. Note that, at least 
in principle, this equation includes possible deviations from the uniformity of 
the asteroid equilibrium temperature on its surface; since the small angular 
sizes of the considered objects, they will be neglected in the following, we will 
work with quantities averaged over the asteroid surface. 
The maximum observed antenna temperature will be obtained when the an- 
tenna beam center points at the asteroid center. By assuming a bi-dimensional 
perfectly symmetric Gaussian beam shape and in the limit of asteroid angular 
sizes much smaller than the beam FWHM, the maximum observed antenna 
temperature can be approximated by the expression: 



rpmax ^ ai n 2T* [R/d] , 

^ iy,ant,obs — u,ant,obj (p^HM/rad) 2 ' 



where T* antobj ~ T eq {e v ~ 1), given by T* antobj = J 47T T U:ant>obj ^,t)dtt/tt obj , 
is the average over the object extent of the object antenna temperature and 
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Qobj = Tc[R/d] 2 is the solid angle subtended by the object, d being its distance 
from the spacecraft and R its effective radius at the considered observation 
time. 

Assuming 150K as mean value of the asteroid antenna temperature, a geo- 
centric distance range between 2 and 4 AU (1.5 x 10 8 km) and a radius range 
between 50 and 500 km, we obtain R/d in the range: 

R/d ~ 2 x 10~ 7 - 2 x 10~ 6 . (8) 

These values can be converted in the maximum antenna temperature, for a 
specific frequency and beam FWHM, using equation (7): 

The antenna temperature range for some PLANCK channels is: 

30GHz , FWHM ~ 33' -> T™^ obs ~ (0.2 - 20)/iK (10) 

100GHz , FWHM ~ 10' -> T™^ obs ~ (2 - 200)^K (11) 

217 - 550GHz , FWHM ~ 5' - T^ obs ~ (10 - 1000);xK . (12) 

Comparing these values with the noise reported in table 2 we can see that, 
mainly at the HFI frequencies, there are good chances to observe asteroids 
and for higher R/d even with the LFI channels. 

In the numerical model we consider only asteroids with R/d larger than 10~ 7 
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(a value referred as "threshold" in the following), since for smaller R/d the 
maximum signal is less than about 1 o- no i S c quoted in Table 2 also for the most 
favorable channels, and there is then no realistic chance to derive information 
on their properties 4 (on the contrary, values of Rj d smaller than this threshold 
have to be considered in the cleaning of Planck data streams from asteroid 
contamination). 

4 Identification of target asteroids in the Planck field of view 

4-1 The numerical algorithm 

To estimate the number of asteroids that will be seen by the LFI and HFI 
channels during the Planck survey and the precise dates of observation we 
have developed a numerical code that calculates the orbits of all large asteroids 
and checks if they fall within the resolution element of a beam. The initial 
orbital elements of all asteroids larger than 50 km are taken from a datafile 
provided by Bowell (2000) while the orbits of the nine planets are derived from 
the JPL ephemerides at the corresponding epoch. The diameters are known 

for most asteroids larger than 50 km and are reported in Bowell's file. For 
4 It can be noted that, since the CMB fluctuations dominate the sky confusion 
noise between ~ 50 and ~ 250 GHz, combining more frequency channels decreases 
the la overall confusion noise quoted for each frequency channel in the last column 
of Table 2. 
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those asteroids whose diameter is not reported in Bowell's file, we used the 
value given in the file of 12,487 asteroids adopted by Zappala et al. (1995) 
for asteroid family search. Whenever the albedo is not directly known from 
observations, an average value, typical of that particular zone of the asteroid 
belt, is used to derive the diameter of the body. To compute the orbits of 
the asteroids and planets we have adopted the RA15 version of the numerical 
integrator RADAU (Everhart 1985) with the highest value of the precision 
parameter. 

At the end of every timestep, we compute the position of each asteroid respect 
to the Earth since the Planck surveyor will orbit together with the Earth 
at its L2 Lagrangian point. Since we know from the mission specifications 
the direction of each beam, we can easily calculate the angle between the 
asteroid and the antisolar direction and the center of the beam of each LFI and 
HFI channel. A possible future refining will be to compute at each epoch the 
position of the spacecraft and compute directly the angle of the asteroid with 
respect to the spacecraft. We fixed the timestep of the numerical integration 
to 1 hr, that is the time required by all the channels to complete a 360° scan. 
Thanks to this choice we can directly compare the position of the asteroid 
with all the channel beam directions during the scan at once. A detection is 
signalled whenever the angle between the asteroid and the antisolar direction 
is within the direction of the beam divided +/— the FWHM/2. 
The input data for the code, apart from the orbital parameters, are the angles 
that define the inclination of each beam respect to the antisolar direction. The 
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code begins the comparison between the position of the asteroids and that of 
the beam at the nominal start of the mission and it continues for two years. 
The code at the end of the run will give the list of those asteroids that will be 
observed by each horn, with the relative date of observation and R/d ratio. 
The output include also a list, for each channel, of the dates when an asteroid 
will be present in the field of view of the beam during the scan. 

4-2 Results 

At the end of the simulation 397 asteroids were observed during the nominal 
length of the mission in at least one of the horns of Planck with a R/d 
high enough for individual detection. For a significant fraction of them (in the 
range from ~ 50 to 100 objects) the good S/N ratio will allow to determine 
the antenna temperature and information on the physical properties of the 
asteroid's regolith. The determination of the antenna temperature at different 
frequencies may be useful to derive models for the surface layers and establish 
some bounds on the depth of the surface layer and the dielectric properties of 
the material. 

The times (i.e. the different consecutive scan circles) an asteroid is consecu- 
tively observed by the same horn during a single passage ranges from 3 to 6 
depending on the FWHM of the channel. It means that on average an aster- 
oid passage lasts about 3-6 hours. In some cases more horns, belonging to the 
same frequency channel, may be crossed within some hours allowing to deter- 
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mine a radio light-curve. In other cases, horns of different frequencies may be 
crossed allowing to get information on the thermal emission of the asteroid. 
From a statistical point of view, asteroids will be observed by each channel of 
both LFI and HFI in 30% of the total scan circles performed in the mission. 
This datum must be taken into account in the extraction of the Galactic and 
extragalactic radio sources from Planck time ordinary data. 
In Fig. 2 the histograms shows the intensity of the signal at detection as mul- 
tiple of the minimum detection threshold R/d ~ 1CT 7 for all the observed 
asteroids on all channels. For large asteroids the radio signal as detected by 
the Planck horns reaches values as large as 12 times the detection thresh- 
old. In the plot at the bottom of Fig. 2 we may notice the signals from the 
5 largest asteroids and, so on at smaller R/d ratio, from all the other visible 
objects. Of course, the number of detected asteroids increases drastically if the 
sensitivity threshold is reduced since the asteroid size distribution is a power 
law. In Table 3 we report the cumulative and differential number of asteroids 
detectable at different levels of the R/d ratio. 

5 Conclusions 

The Planck ESA mission will perform a high-angular resolution mapping 
of the microwave sky over a wide range of frequencies, from 30 to 900 GHz. 
During the two years planned for the mission the surveyor will frame a large 
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Fig. 1. Number of observations of asteroids by Planck, during the whole mission at 
different sensitivity levels (xlCT 7 ). In the bottom plot at high signal intensity only 
the 5 largest asteroids appear in Planck's horns. The increasing number of small 
asteroids that can be observed at lower sensitivity explains the growing amount of 
detections. 
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number of main-belt asteroids. In this paper we investigate in detail for the 
first time the possibility for Planck to detect main-belt asteroids. The main 
parameter in the asteroid detection process by the Planck horns is the R/d 
ratio. The instrumental sensitivity to such objects has then been denned as 
the minimum R/d value required to safely observe them. This value is given 
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Table 3 



Cumulative and differential number of asteroids may be detected at different R/d 
values. 



R/d 


Differential number 


Cumulative number 


1-2 xl(T 7 


299 


397 


2-3 xl(T 7 


76 


98 


3-4 xl(T 7 


15 


22 


4-5 xl(T 7 
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5-6 xl(T 7 





3 


6-7 xl(T 7 





3 


7-8 xl(T 7 





3 


8-9 xl(T 7 


1 


3 


9-10 xl(T 7 


1 


2 


1-2 xl(T 6 





1 


2-3 xl(T 6 





1 


3-4 xl(T 6 


1 


1 



by the Planck sensitivity at different frequencies rescaled to take into ac- 
count the integration time and the confusion noise due to the background. 
An accurate evaluation of the noises and the sensitivity degradation, related 
to the high proper motion of the asteroids, yielded a minimum R/d ratio of 

(R/d) mhi ~ icr 7 . 

A numerical simulation of a two-years mission has been performed in order to 
estimate the number of objects whose (R/d) is greater than about 1CT 7 and 
that will then be observed by Planck. The simulation uses updated cata- 
logs for the orbital elements of the main-belt asteroids and their diameters. 
The orbit of each asteroid has been integrated with a very short timestep to 
compute the relative position of the body with respect to the Planck horns 
and the value of (R/d). An accurate mission simulation has been obtained. 
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Up to 397 asteroids are expected to be detected in the various Planck chan- 
nels (Table 3). Detectable asteroids will appear in about 30% of the total sky 
circles scanned by the mission. 

The previous survey of Main Belt Asteroids was performed by the Infrared 
Astronomical Satellite (IRAS) on 1983 in four wavelength bands centered near 
12, 25, 60 and 100/im, much lower than the minimum wavelength observed 
by Planck of 350/zm, corresponding to the highest frequency of 857 GHz. It 
surveyed approximately 96% of the sky and 2228 different multiply observed 
asteroids were associated to IRAS sources (Tedesco et al. 2002), providing a 
good estimate of diameter and albedo for most of them. 
IRAS's 12/iin limiting sensitivity, for S/N=3, was about 150 mJy (Tedesco 
and Desert 2002), and assuming a Main Belt Asteroid as a black body at 150 
K this flux can be translated to about 56 mJy at 350/xm (857 GHz), well below 
the Planck noises reported in Table 2. 

PLANCK will provide flux measurements for a smaller sample of asteroids 
(~ 50 — 100 objects) compared to IRAS, but at different wavelengths, al- 
most unexplored for this class of Solar System objects. This will improve our 
understanding of the thermal emission and the related surface properties of 
asteroids. 
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